Proceedings Title lAU Symposium 
Proceedings lA U Symposium No. 252, 2008 
Licai Deng, K.L. Chan & C. Chiosi, eds. 



© 2008 International Astronomical Union 
DOI: 00.0000/XOOOOOOOOOOOOOOOX 



Around the Pair Instability Valley - Massive 

SN Progenitors 

Roni Waldman 

Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 
email: waldman@cc.huji.ac.il 



oo 

o 

O 

>>■ 



Oh 
I 

O 



> 

00 

rn 

o 

oo 

o 



- 

X 



Abstract. The discovery of the extremely luminous supernova SN 2006gy, possibly interpreted 
as a pair instability supernova, renewed the interest in very massive stars. We explore the 
evolution of these objects, which end their life as pair instability supernovae or as core collapse 
supernovae with relatively massive iron cores, up to about 3Mq. 
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1. Introduction 

The interest in the evolution of very massive stars (VMS), with masses > 100 Mq, 
has recently b een revived by the discovery of SN 2006gy - the most luminous supernova 
ever recorded ( Smith et al. l2007l ). This object, having a luminosity of ~ 10 times that 
of a typical core-collapse SN (CC SN), is probably the first evidence of a pair instability 
SN (PISN) IWooslev et aP (|20n7l ). PISN are massive stellar objects, whose evolutionary 
path brings their center into a region in thermodynamical phase space (p < 10^, T > 
10^), where thermal energy is converted into the production of electron-positron pairs, 
thus resulting in loss of pressure an d hydrodynamic inst ability. This type of supernova 
was first suggested 4 years ago by iBarkat et al.l (119671). and since then several works 
were carried out (e.g. lFralevlll968t fober et al.lll983tlEl Eid et al.lll983HBond et al.lll984 : 



Heger fc WooslevI 120021: Iffirschi et al.ll2004HEldridge fc Toutll2004 iNomoto et al.H2005[ ). 
however the overall interest in this topic has been relatively small, mainly due to lack of 
observational data. 

It was originally believed that stars massive enough to produce PISN could only be 
found among population III stars with close to zer o metallicity (Z 5, 1 0^"^), and hence 
only at very high redshift (z > 15). More recently IScannapieco et aL (20ol) discussed 
the detectability of PISN at redshift of z ^ 6, arguing that metal enrichment is a local 
process, ther e fore m etal-free star- forming pockets may be found at such low rcdshifts. 
Langer et al.l ( 20071 ) introduced the effect of rotation into studying this question conclud- 
ing that PISN could be produced by slow rotators of metalli city Z < Z(^/3 a.t a rate of 
one in every 1000 SN in the local universe. Furthermore, Smith et al. ( 2007? ) point out, 
that mass loss rates in the local universe might be much lower than previously thought, 
so that massive stars migh t be left w i th eno ugh mass to become PISN. This conclusion 
is also supported by Yungelson et al.l (|2008l ) who extensively discuss the mass loss rates 
and fates of VMS. It is interesting to note, th at SN 2006gy took place in the nearby 
Universe. Following the discovery of SN 2006gv. 'Um eda fc Nomotol ( 2008) addressed the 
questi on of how much ^^Ni can be produced in massive CCSN, while Heger fc WooslevI 
( 2008[ ) computed the detailed nucleosynthesis in these SNe. 

The interest in VMS is further motivated by the discovery of Ultraluminous X-ray 
Sources (ULX), which can be interpreted as mass-accreting intermediate mass black 
holes (IMBH) with mass - (10^ - 10^) Mq. One of the possible scenarios for IMBH 
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formation is by VMS formed by stellar mergers in compact glo bular clusters (see e.g . 
Yungelson et al. 20081 and references therein). In this context, Nakazato et al.l ( 20061 



2007? ) studied the collapse of massive iron cores with M >3 Mq. In their first paper they 



treat the fate of stars of mass ^ 300 Mq which reach the photodisintegration temperature 
(w 6 X 10^i<') after undergoing pair instability. The entropy per baryon of these models 
at photodisintegration is s > IGfcs compared with the classical core-collapse SN with 
s ~ Iks- In the second paper they aim to bridge this entropy gap, corresponding to core 
masses of (3 — 30) Mq but claim that there is a lack of systematic progenitor models for 
this range, hence they use synthetic initial models for their calculations. 

In this work we focus mostly on the mass range M < 80 AIq (He core mass Mhc ^ 
36 Mq) immediately below the range which enters the pair instability region, and present 
a systematic picture of the resulting CCSN progenitors. 



2. Method 

Since the mass loss rat es of stars in this range are highly uncertain, (see e.g. discussion 



b vlYungelson et al 



20081). we avoid dealing with this question by following the example 



of lHeger fc Wooslev ( 2002 ). and modeling the evolution of helium cores. Our helium core 
initial models are homogeneous polytropes composed entirely of helium and metals, with 
metallicity Z « 0.015, in the mass range (8 — 160) M© . The models were then evolved 
to the helium zero age main sequence. In the following we will refer to these models as 
"HeAf" where N is the mass of the model. For comparison we evolved also a few models 
of regular hydrogen stars, beginning from the zero-age main sequence (ZAMS). We will 
refer to these models as "MA^" where N is the mass of the model. All our models have no 
mass-loss. We argue that as long as the mass loss rate is not so high that it will cut into 
the He-core, the evolution after the main-sequence phase will be virtually independent of 
the fate of the hydrogen-rich envelope. We followed the evolution of each model until the 
star is either completely disrupted (for the PISN case) or Fe begins to photo-disintegrate 
(for the CCSN case). 

We followed the evolution using the Lagrangian one dimensional Tycho evolutionary 
code version 6.92 (with some modi fications), publicly available on the web (the code is 
described in lYoung fc Arnett 2005 ). Convection is treated using the well known mixing 



length theory (MLT) with the Ledoux criterion. In the MLT formulati on of Tycho, the 
value of the mixing length parameter fit to the Sun is umlt ~ 2.1 ( Young fc Arnet"3 



2005h . so we used a value of umlt = 2 in our calculations. The nuclear reaction 



rates used by T YCHO are ta ken from the NON-SMOKER database as described in 
iRauscher fc ThielemannI (|2000l ). 



3. Results and Discussion 

Among the He-core models we computed, those in the mass range Al ^ 36 Mq do 
not reach pair instability and end their lives as CCSN. Fig. [T] shows the evolution of 
the central density and temperature (left panel) and the density structure of the pre-SN 
(right panel), at the moment when the central temperature reaches 7 x 10^ if. The two 
extreme models He8 and He36 are shown, as well as M80 which has a He-core mass 
similar to the He36 model, and M20 - a typical CCSN progenitor. The left panel also 
shows two He-core models - He80 and Hel60 that reach pair instability. Note that the two 
models He36 and M80 are indeed very close to each other. An example of the composition 
at the pre-SN stage (Tc = 7 x IQ'^K) is shown in Fig. O The Fe-core mass is « 3Mo, 
topped by a shell of k, 10 Mq of Si-group elements. The size of the Fe-core (defined as 
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Figure 1. Evolution of tlie central density and temperature (left panel) and pre-SN density 
structure (right panel). Each line is labeled "M" for stellar models and "He" for He-core models, 
followed by the mass of the model. 

the mass coordinate where the electron mole fraction Ye < 0.49) and the central entropy 
per baryon at the pre-SN stage for the whole set of models is shown in Fig. [3l Note that 
the size of the Fe-core is slightly non-monotonic. The central entropy, is monotonic with 
mass, but slightly differs between He-core and stellar models. 

As can be seen in the above figures, the outstanding features of these massive models 
compared with the typical CCSN example M20 are: 

(a) Lower central density and higher central entropy. 

(&) A much shallower density profile. 

(c) Relatively large Fe-cores, up to about SAf©, and a large amount (up to about 
10 Mq) of Si-group elements. 

These differences might have a considerable impact on the behavior of these models 
during core collapse and on the outcome of the explosion, a question which we hope to 
address in the future. 
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Figure 2. Pre-SN composition of tlie He36 
model. "Si" and "Fe" stand for tlie total of Si- 
and Fe-group elements respectively. 



Figure 3. Mass of the Fe-core (squares) and 
central entropy per baryon (circles) for the 
computed models. Filled shapes designate He- 
core models, open shapes - stellar models. 
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